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Schwinger's Dynamical Casimir Effect is one of several candidate explanations for sonoluminescence. 
Recently, several papers have claimed that Schwinger's estimate of the Casimir energy involved is 
grossly inaccurate. In this letter, we show that these calculations omit the crucial volume term. 
When the missing term is correctly included one finds full agreement with Schwinger's result for 
the Dynamical Casimir Effect. We have nothing new to say about sonoluminescence itself except 
to affirm that the Casimir effect is energetically adequate as a candidate explanation. 
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I. INTRODUCTION 

Several years ago, Schwinger wrote a series of pa- 
pers wherein he calculated the Casimir energy re- 
leased in the collapse of a spherically symmetric bubble 
or cavity. Using general arguments, he showed that the 
effect was mostly a volume effect, and derived a simple 
and elegant formula for the energy release involved in 
the collapse. He found that (for each polarization state) 
the "dielectric energy, relative to the zero energy of the 
vacuum, [is given] by 



E = -V 



(1) 



5*0 the Casimir energy of a uniform dielectric is nega- 
tive" . From the above one finds that a dielectric slab 
with a spherical vacuum cavity of radius R has a higher 
Casimir energy than the same slab of material with 
the cavity re-filled with dielectric. Introduciag a wave- 
number cutoff K into the previous expression, □ and sum- 
ming over polarization states, shows that the Casimir en- 
ergy of a cavity in a dielectric, relative to pure dielectric, 
is 
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In general, this volume term will be the dominant con- 
tribution. 



In view of the elegance and simplicity of this result, it 
is natural to ask whether it can also be derived by more 
traditional quantum field theoretic means. Indeed, the 
existence of such a volume contribution is easy to verify 
on general physical grounds: 

(1) One can view Schwinger's result in elementary terms 
as simply the difference in zero-point-energies, obtained 
by integrating the difference in photon dispersion rela- 
tions over the density of statesLl 
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At low wave-numbers, we know that the dispersion re- 
lation for a dielectric is simply summarized by the zero- 
frequency refractive index n. That is 



LL!{k) — > ck/n 



as 



fc 0. 



(4) 



On the other hand, at high enough wave-numbers, the 
photons propagate freely through the dielectric: They are 
then simply free photons travelling through the empty 
vacuum between individual atoms. Thus 



uj(k) — > ck 



as 



k 



(5) 



From the above we know that the integrand must go to 
zero at large wave-number. In fact for any real dielectric 
the integrand must go to zero sufficiently rapidly to make 
the integral converge, since after all we are talking about 
a real physical difference in energies. 

To actually calculate this energy difference one requires 
a suitable physical model for uj{k). Schwinger's calcula- 
tion is equivalent to picking the particularly simple 
model 



*For sonoluminescence, this wave-number cutoff can be re- 
lated to the wavelength of the electromagnetic radiation emit- 
ted in the collapse of the bubble. For generic dielectrics, 
this wave-number cutoff is related to the high-wave-number 
asymptotic behaviour of the dispersion relation. 



^The dots denote finite- volume corrections. We shall develop 
this density-of-states point of view more fully in a separate 
publication. 
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uj{k) ^ — e{K -k) + ck e(fc - K). 



(6) 



II. PHYSICAL DESCRIPTION OF THE 
CALCULATION 



Here <d{x) is the Heaviside step function, and K is 
a wave-number which characterizes the transition from 
dielectric-hke behaviour to vacuum-hke behaviour. Note 
that the cutoff K describes an actual physical situation: 
It is a surrogate for all of the complicated physics that 
would be required to make a detailed model for the di- 
electric to vacuum transition. 

(2) We also know that the quantum action in 3+1 dimen- 
sions generically contains divergences which range from 
quartic to logarithmic, in addition to finite contributions. 
As is well known, this "cosmological constant" contribu- 
tion (the quartic divergence) will not vanish unless the 
theory has very special symmetries (for example — super- 
symmetry). Thus energy densities that go as (cutoff)^ 
are generic in (3+1) dimensions. 

(3) Alternatively, one could perform an explicit quan- 
tum field theoretic calculation of the Casimir energy in 
some model problem and thereby verify Schwinger's re- 
sult. A step in this direction has been provided by Mil- 
ton et al. ^,||,D , who attempted to compute the Casimir 
energy associated with a spherical cavity of radius i?, di- 
electric constant ei, and permeability //i, embedded in 
an infinite medium with dielectric constant £2 and per- 
meability ^2- They found that the dominant term is not 
the volume term but a surface term which is proportional 
to R^K^{ei — £2)^- We, however, have re-analysed these 
calculations and do find a volume term which dominates 
except for very small bubbles. 

We have performed the calculation of the Casimir en- 
ergy in two different and complementary ways: 

(a) We have taken the formalism of Milton et al. |^,||j6| 
and applied it directly to an ab initio calculation of the 
Casimir energy. We compute the energy difference be- 
tween the following configurations: (Case I) an otherwise 
uniform medium with dielectric constant £2 and perme- 
ability /i2 , containing a spherical cavity of radius R with 
dielectric constant £1 and permeability /ii, and (Case II) 
a completely uniform medium with dielectric constant £2 
and permeability fi2- This energy difference is given as 
a sum over a series of integrals involving Ricatti-Bessel 
functions. Some of the sums can be evaluated explicitly 
while others can only be evaluated by using an asymp- 
totic analysis of the type used by Milton et al. We verify 
the existence of both volume and sub-dominant surface 
contributions. 

(b) We have analyzed the extant calculations to see ex- 
actly where they differ from the present calculation. We 
find that the subtraction scheme they use to calculate 
the Casimir energy does not correspond to the physical 
situation in question. We isolate the difference in en- 
ergy between these calculations and the correct one. We 
will explicitly show that this difference is proportional to 
volume. 



Milton et al. explicitly calculated the electro- 

magnetic Green functions for a dielectric ball embedded 
in an infinite space of (different) dielectric material, and 
then attempted to calculate the Casimir energy by explic- 
itly integrating these Green functions over "all space" . 
Note that an important limitation of any such calcula- 
tion is that any attempt at explicitly calculating Green 
functions must be restricted to systems of extremely high 
symmetry — such as half-spaces, slabs, or balls. The basic 
strategy is to take the classical expression for the energy 
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promote the electric and magnetic fields to be operator 
quantities, and then calculate the vacuum expectation 
value 

E=\l [e{E{Q,x)-E{0,x)) 

^ J Geometry 

+ ^^{B{0,x)-B{0,x))]£x. (8) 

The geometry is incorporated in the calculation both via 
the limits of integration and via the boundary conditions 
satisfied by the fields. Since these two-point functions 
are of course divergent, they must be rendered finite by 
some regularization prescription. Milton et al. use point- 
splitting in the time direction: 

E{t)^\[ [e{E{T,x)-E{0,x)) 

^ J Geometry 

+ j^{B{T,x)-BiO,x))]d^x. (9) 

All of the technical aspects of the analysis then focus 
on the calculation of these two-point correlation func- 
tions (Green functions) by explicitly solving for the TE 
and TM modes appropriate for a spherical ball with di- 
electric boundary conditions; and then explicitly writing 
down the Green functions as a sum over suitable combi- 
nations of Ricatti-Bessel functions and vector spherical 
harmonics. To avoid unnecessary notational complica- 
tions, we schematically rewrite the above as 

£;(r) = i / G[,,^](T,x;0,x)d3x, (10) 

^ J Gco7netry 

where G[g^^] {t, x; t' , x') is simply shorthand for the linear 
combination of Green functions appearing above. 

We may calculate these Green functions for three dif- 
ferent geometrieaa: 



■'■Notice that in Milton et al. the dielectric properties of these 
media are taken to be frequency independent, the cutoff being 
put in "by hand" via time-splitting. 
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Case I: A dielectric ball of dielectric constant ei, perme- 
ability /zi , and radius R embedded in a infinite dielectric 
of different dielectric constant 62 and permeability fi2- 
(In applications to sonoluminescence, think of this as an 
air bubble of radius R in water.) 

Case II: A completely homogeneous space completely 
filled with dielectric (e2,/i2). (In applications to sonolu- 
minescence, think of this as pure water.) 
Case III: A completely homogeneous space completely 
filled with dielectric (ei,/zi). (In applications to sonolu- 
minescence, think of this as pure air.) 

We are in complete agreement with the extant cal- 
culations and results for these three individual Green 
functions — where we disagree, as will be shown, is in the 
way that these three Green functions are inserted into the 
computation for the Casimir energy. 

Milton et al. calculate an "energy difference", which 
we will call Egurfacei and which they define as 



E sur f ace 



Gr 



all : 



G 



r>R 



GlII 



r<R 



(11) 



The computation of this quantity in [^H^,^ is mathemat- 
ically correct. An asymptotic analysis shows that this 
expression is indeed proportional to the surface area — 
plus even higher-order terms. However, the energy cal- 
culated from the above expression does not correspond 
to the energy of a physically realizable situation. 
The physically correct quantity to compute is ||l| 
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(12) 



Observe that this quantity is simply the difference in 
energy between two real physical situations: (Case I) 
having the dielectric ball present and (Case II) replac- 
ing the dielectric ball by the surrounding medium. This 
is exactly the quantity that Schwinger calculates in ref- 
erence JJ/ to describe the Casimir enegy released in the 
collapse of the bubble: it is the energy released in evolv- 
ing from bubble to no-bubble. The difference between the 
two calculations is 



AE = E, 



Casimir ^surface 



o / {Giii - Gji} . 

^ Jr<R 



(13) 



This difference is easily seen to be proportional to the 
volume: remember that G/77 and G// are Green func- 
tions corresponding to two spaces that are completely 
filled with homogeneous dielectrics — therefore they are 
each individually translation invariant. (When one ex- 
presses these Green functions in terms of spherical polar 
coordinates this is not obvious.) This observation per- 
mits one to pull the Green functions outside the integral, 
so that 



AE = -V {Gin{T, 0; 0, 0) - G//(t, 0; 0, 0)} , (14) 

where V is the volume of the ball of radius R. We shall 
now show that this term is in fact exactly in conformity 
with Schwinger's result. 



III. "AB INITIO" CALCULATION 
A. The energy density 

We now calculate the Casimir energy for the geomet- 
rical configuration previously described. We use tech- 
niques developed by Milton et al., but use, as Schwinger 
did, a wave number cutoff and shall present the calcula- 
tion in as much detail as space permits. We defer many 
technical details to a forthcoming publication. 

For each individual geometry, the energy density T" 
can be evaluated by using the dyadic Green function for- 
malism (Henceforth we use natural units.) For 
Case I one finds 



Tf (r) = Re 



(15) 



with identical expressions holding for the other geome- 
tries. Here the w-integral arises from the time-splitting 
regularization, while we have used the notation k = \uj\n 
with n the appropriate position-dependent refractive in- 
dex (ni inside the dielectric sphere, 712 outside), and 
have defined the quantity A/(fc,r) by 



Xi{k,r) = Y,{2i+\)< 



2 , ^(^+1) 



F/(/c;r, r) 
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^r^d^.^'d^/' [^/(fc;-i>-^)] 

+ [{fD ^ {G[)] . 



(16) 



The functions Fiir^r') and Ge{r,r') are the Green func- 
tions for the electrical and magnetic fields in the ap- 
propriate geometry, and are given below. Similar re- 
sults, with appropriate substitutions for the momenta, 
hold when one makes reference to Cases II and HI. 
Note that when making the substitutions (/) {II) or 
(/) (III), one should also change the refractive index 
that implicitly appears in the factor k. In addition it 
should be borne in mind that fc/ is a function of position: 
kj ~ ni\Lu\ ~ km inside the dielectric sphere, whereas 
kj — n2\Lo\ — kji outside the dielectric sphere. 

Because Fi and Gi depend only on the absolute value 
of Lo we can write the energy density as 



Tf (r) = Re 
= Re 
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(17) 



The Casimir energy, Eq. ( [12D is obtained by taking the 
difference in energy densities and integrating over aU 
space while paying attention to the appropriate index 
of refraction for each region of space. Thus 



Casimir 
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xIni[X/(fc,r)-X//(fc,r)]. 



(18) 



This expression for the Casimir energy is completely 
equivalent to equation (41) of ||], and equation (4.2b) 
of 1^ and is also closely related to equations (30a) and 
(30b) of Q. Note that extant calculations use the same 
time-splitting parameter for the two different media — the 
physics behind this choice is far from clear, and we shall 
return to this point in a future publication. 

It is now clear how one should modify these expressions 
to replace time-splitting regularization by a wave-number 
cutoff. For generality we can take an arbitrary wave- 
number cutoff described by some smooth real function 
f{k) which goes to zero as fc — + oo and simply write 
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xlm[f{ki)Xi{k, r) - f{kii)Xii{k, r)]. (19) 



With due caution, the relevant Green functions can be 
read off from 



B. The Green functions 



The function jt [x) is the spherical Bessel function of or- 
der ^ and hi{x) = h'^p {x) is the spherical Hankel function 
of the first kind. See equations (12a) and (12b) of [Q, 
equation (16) of P, or equation (2.13) ofJOl. The quan- 
tities A^p Q and q are those given in . 
Case II: 
For all ri, 

FP, G'/iri,r2) = ikii ]i{kiir<) h,{kiir>). (24) 

Case III: 

For all ri, 

F/",Gf ^(ri,r2) = ikui jdkiur<) hi{kuir>). (25) 

We are now ready to explicitly compute the Casimir en- 
ergy. In passing we remark that the object F^^^ defined 
in equation (29) of Q|, equation (35) of and equation 
(3.7) of j^, which is essential for those calculations, is not 
a Green function of any differential operator. Specifically, 
Ff'^ does not satisfy the dielectric boundary conditions. 
It is not even continuous, and is merely a potpourri of 
two different Green functions which does not have any 
particular physical relevance. 



C. The Casimir energy 

We calculate Ecasimir using equation (|l^), [equiva- 
lently (|l|)] together with equations (|2|-|2|). When eval- 
uating the imaginary parts of X it is convenient to in- 
troduce the Ricatti-Bessel functions Si{x) = xj^{x) and 
ee{x) = xhi{x). One also needs the identity 



In evaluating the Green functions one must be careful 
to correctly incorporate the boundary conditions appro- 
priate to the geometry and the physics. This means that 
they must satisfy appropriate continuity conditions de- 
rived from Maxwell's equations. That is 



-Si{x) = s'l{x) + Si{x), 



E±, eEr, and Br 



(20) 



must be continuous. In terms of and G^, one sees that 

d 1 ^ 

fiFi, Ge, —rFe, and -^rGe, (21) 
or e or 

must be continuous. The Green functions are 
Case I: 

For ri, r2 < R: 

Ft. iG\{ri,r2) = ikm je{kHir<) 

X [he{kiiir>) - A^^ jiikmr^)] . (22) 

For ri,r2 > R'- 

F/,G|(ri,r2) = iku ht{kiiry) 

X []i{kiir^) - B'pc heikiir^)] . (23) 



(26) 



together with an identical equation which holds for ee{x). 
After some rearrangement (technical details are sup- 
pressed and will be relegated to a more detailed forthcom- 
ing publication) we find that inside the dielectric sphere 

Im{X/(fc,r)}i„ 

i=i 

X {2[si{x)]^ + [4{x)]^ + Se{x)s';{x)} \k,,,r 

-tlRY.i^i+l)ne{A'^ + A%} 

^ e=i 

X {2[seix)f + [s',ix)f + seix)s';ix)} 

(27) 

A remarkable Ricatti-Bessel function identity permits 
us to perform the first sum over i exactly. Using 

oo 

^(2^ + 1) [2si{xf + s'i,{x)^ + se{x)s';{x)] = 2x^ , (28) 
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one obtains that 

Im{X/(fc,r)}i„ = 4k^jjj 
km 



1=1 

y.{2[s,{x)f + [s',{x)f + s,{x)4{x)]l^^^^. 

(29) 

Note that for the Q-terms we cannot exphcitly per- 
form the £ summation because of the comphcated I de- 
pendence of these coefficients 

It is very important to notice at this point that the 
A^p Q-terms in the above expressions are the only pieces 
retained in the currently extant calculations, the other 
terms unfortunately have been missed there due to the 
use of the wrong "Green functions" . 

Taking a cue from the above, the results for the region 
inside the bubble can be written (using self explanatory 
notation) as 

lni{Xi{k, r)}„ = 4n?|c.|3 + Q,^'->-(fc„,, r). (30) 

Turning to the region outside the dielectric sphere, one 
gets that 

luv{Xi{k,r)}out 

J OO 

X {2[si{x)f + [s'f(x)f + se{x)s';{x)} \k„r 
e=i 

x{2[ee{x)r + [e',{x)r+e,{x)enx)}]\,^^^^ 

J OC 

kjj 



4klj 



^5:(2^+l)Re[(B|^ + i?S) 

e=i 

x{2Mx)r + [e',{x)r+e,{x)e';{x)}]l^^^^ 
Again, in self-explanatory notation 

lm{Xiik,r)Ut = 4n3|c.|3 + Q— /-(fc,,,r). 
For Case II one simply has 



(31) 



(32) 



X [nlf{ni\uj\) - nlf{n2\uj\)] 
+ I ^'dr I |^/(ni|c.|)g->-(fc,r) 



^0 

OO 7 

2- / -^/(n2|c.|)g->-(fc,r). 

R Jo 



-f / r^dr 



(34) 

The remaining integrals for the n^|ci;|^ term are triv- 
ial. Changing the integration variable to fc = and 
explicitly re-inserting the appropriate factors of Ti and c, 
we get 
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which is the central result of this paper. 
This is explicitly of the form: 

(Schwinger's volume term) + (surface term). 

The "surface term" corresponds to Esurface and is given 
by the two double integrals in the expression for Ecasimir 
above. For time-splitting regularization and dilute dielec- 
tric media, these terms were explicitly shown by Milton 
et al. to be proportional to the surface area (plus even 
higher-order corrections). The first term is the volume 
term not present in some of the existing calculations. In 
fact, after approximating air by vacuum (setting ni = 1) 
and using Schwinger's momentum space cutoff, this inte- 
gral is exactly equal to Schwinger's result [Q. 



IV. DISCUSSION 

The main result of this paper can be succinctly stated: 
in a dielectric medium of dielectric constant n the 
Casimir energy of a cavity — the difference in zero point 
energies of a dielectric medium of refractive index n with 
and without a vacuum cavity of volume V — is: 



lm{Xii{k,r)}aU space = 2 ^(2£+ 1 



e=i 



X {2[se{x)]^ + [4(x)]2 + se{x)s';{x)} \k,,r 



(33) 



Going back to the momentum-space regulated Casimir 
energy, equation (|l^), we obtain 
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r^dr 



dhJ 

2^' 
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Ecasimir = V HcK'^ 



1 



1 



(36) 



with this volume term dominant if the scale of the bubble 
is larger than the cutoff wavelength 2tt/K. This result 
is completely in agreement with Schwinger's calculation 
in [0, and Schwinger's argument is now buttressed by 
our explicit re-assessment of Milton et alJs calculation 
for a spherical dielectric ball. 

We close with what is perhaps a minor point that we 
nevertheless feel should be made explicit: the volume 
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contribution to the Casimir energy is always there, and 
is always physical, but it is sometimes safe to neglect it. 

For example, a situation equally physical as the one 
we have considered here is the following: suppose one is 
provided with a fixed number of dielectric bodies of fixed 
shape (in particular, of fixed volume), and suppose that 
one simply wishes to move the bodies around in space 
with respect to each other. Then the bulk volume con- 
tributions to the Casimir energy, while still present, are 
constants independent of the relative physical location of 
the dielectric bodies, and so merely provide a constant 
offset to the total Casimir energy. // all we are interested 
in is the energy differences between different spatial con- 
figurations of the same bodies then the various volume 
contributions can be quietly neglected. 

On the other hand, the volume contribution is of crit- 
ical importance whenever one wants to calculate the en- 
ergy difference between an inhomogeneous dielectric and 
a homogeneous dielectric wherein the irregularities have 
been filled in. This is, precisely, the physical situation in 
the case of bubble collapse in a dielectric medium. 
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